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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  results  of  the  research  were 
obtained  by  ARO,  Inc.  (a  Sverdrup  Corporation  Company),  AEDC  Division,  operating 
contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under  ARO  Project 
Number  Y32K-15.  Mr.  Marshall  Kingery,  AEDC/DOT,  was  the  Air  Force  project  manager. 
The  manuscript  was  submitted  for  publication  on  September  24,  1979. 
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1.0  INTRODUCTION 

There  is  a  continuing  need  to  analyze  the  molecular  flow  in  space  chambers,  through 
baffles,  in  ducting,  and  in  similar  complex  systems.  When  chemical  or  ion  rockets  are  tested 
under  high  altitude  space  conditions,  it  is  desirable  to  minimize  the  pressure  or  backflow  in 
certain  regions  of  the  test  chamber.  Because  of  the  high  velocity  of  molecules  emerging  from 
a  rocket  exhaust  and  the  fact  that  these  molecules  are  not  immediately  accommodated  to 
chamber  conditions,  a  systematic,  analytic  approach  to  the  pretest  chamber  configuration 
optimization  is  desirable. 

The  General  Unwanted  Energy  Rejection  Analysis  Program  (GUERAP)  has  been  used 
extensively  in  the  past  few  years  to  analyze  a  variety  of  optical  systems.  GUERAP  (Ref.  1) 
uses  a  Monte  Carlo  ray  trace  technique  to  analyze  the  radiation  characteristics  of  complex 
optical  systems. 

There  are  many  similarities  between  radiation  ray  trace  and  molecular  ray  trace 
techniques.  It  was  decided  that  the  most  expeditious  approach  to  providing  the  molecular 
flow-field  analysis  would  be  to  utilize  the  geometric  and  ray  trace  portion  of  GUERAP.  This 
modified  GUERAP  program  is  called  the  General  Unwanted  Molecule  Rejection  Analysis 
Program  (GUMRAP). 

The  major  modification  to  GUERAP  was  the  elimination  of  several  unnecessary 
radiation  subroutines  and  the  addition  of  a  gas-surface  interaction  model.  Shown  below  are 
the  changes  made  in  GUERAP. 

1.  Removal  of  diffraction  of  rays 

2.  Removal  of  refraction  regions 

3.  Removal  of  the  important  surface  technique 

4.  Removal  of  internal  radiation 

5.  Removal  of  irrelevant  radiation  scattering 

6.  Addition  of  realistic  gas-surface  interactions  model 

Time  limitations  precluded  the  elimination  of  all  FORTRAN  statements  and  variables 
that  are  remnants  from  GUERAP. 
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2.0  GAS-SURFACE  INTERACTION  MODEL 


It  has  been  100  years  since  Maxwell  proposed  that  when  a  gas  reflects  from  a  solid 
surface  a  certain  fraction  of  the  molecules  is  reflected  diffusely  (fully  accommodated)  and 
the  remainder  of  the  molecules  are  reflected  specularly  (no  accommodation).  Over  the  years 
other  investigators  such  as  Nocilla  (Ref.  2)  and  Epstein  (Ref.  3)  have  proposed  more  realistic 
gas-surface  interaction  models.  However,  all  previous  models,  while  being  more  or  less 
successful  in  explaining  energy  exchange  with  a  surface,  failed  to  realistically  describe  the 
scattering  of  individual  molecules. 

More  recently,  Cercignani  (Ref.  4),  Kuscer  et  al.  (Ref.  5),  and  Wenaas  (Ref.  6)  have 
shown  that  any  gas-surface  interaction  model,  in  addition  to  satisfying  the  obvious 
condition  of  maintaining  the  equilibrium  velocity  distribution  at  the  surface  temperature, 
must  satisfy  the  reciprocity  relation.  The  introduction  of  a  scattering  kernel,  P(v',v),  which 
gives  the  reflected  probability  density  as  a  function  only  of  incident  and  reflected  velocity 
vectors,  implicitly  assumes  that  a  molecule  interacts  with  the  wall  on  an  individual  basis, 
independent  of  other  molecules.  In  other  words,  a  molecule  with  velocity  “v”  does  not 
“know”  whether  it  is  a  rocket  exhaust  or  part  of  a  gas  in  equilibrium  with  the  wall.  This 
result  permits  the  reciprocity  condition  on  the  scattering  kernel,  which  was  derived  for 
equilibrium  conditions,  to  be  applied  for  nonequilibrium  conditions. 

2.1  BASIC  EQUATIONS 

The  time-independent  interactions  between  a  gas  and  a  surface  may  be  expressed 
mathematically  by  means  of  a  scattering  kernel,  P(v',v),  which  relates  the  probability  that  a 
molecule  with  velocity  v'  will  be  scattered  by  a  surface  into  a  velocity  close  to  v  (refer  to 
Fig.  1). 

Kuscer  et  al.  (Ref.  5)  point  out  the  importance  of  reciprocity  in  connection  with  a 
mathematical  description  of  gas-surface  interaction.  The  reciprocity  relation  which  must  be 
satisfied  by  any  scattering  kernel  is 

—  Vy  P(v  v)  exp  (—  v'V  2RTw)  =  VyP(—  v,—  v  0  e  xp  ( — v2/2RTw)  (1) 

In  addition,  P(v',  v)  must  satisfy  the  obvious  condition  of  being  non-negative  and,  for  no 
mass  exchange  with  the  surface,  normalized  in  half  velocity  space.  In  other  words, 
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and 
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/ 


oo 


/ 


oo 

f  P  (vVv)  dv  dv  dv  =  1 

J  X  )  z 

—  oo 


(3) 


Equation  (3)  states  the  conversation  property  of  molecules,  i.e.,  that  molecules  are  assumed 
not  to  be  trapped  or  cryopumped  by  a  surface. 


y,vy 


Figure  1.  Coordinate  systems  for  gas-surface  interaction. 

Kinslow  (Ref.  7)  has  developed  a  product  expression  for  P(v',  v)  that  satisfies  the 
conditions  expressed  in  Eqs.  (1),  (2),  and  (3)  and  also  contains  the  experimentally  observed 
characteristics  of  gas-surface  interactions  as  discussed  in  Refs.  7,  8,  and  9. 


The  scattering  kernel  presented  in  Ref.  7  is  as  follows: 

2 A2  A  A 


P  ft  'ft) 


y  x  z 


(2RT 


S  I  (-2 A  VA2-1  S'S  ) 
3/  y  0  y  y  v  r 


-p  s;-A*sJ-(A2-i)s;2-A2sy2-(VF7T  s;-azs^!] 


(4) 
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where  Ax,  Ay,  and  Az  are  arbitrary  constants  greater  than  1  and  are  related  to  the  gas-surface 
energy  accommodation  coefficient,  a,  by 


A-2 


=  ay,  etc. 


(5) 


For  the  special  condition  of  ax  =  ay  =  az  =  1  (fully  accommodated),  Eq.  (4)  reduces  to 
a  fully  accommodated  Maxwellian  velocity  distribution  at  the  surface  temperature,  with  a 
cosine  spatial  reflection.  For  the  condition  of  ax  =  «y  =  az  =  0  (no  accommodation),  Eq. 
(5)  gives  specular  reflection.  These  correct  limits  help  support  the  validity  of  the  scattering 
kernel  as  given  by  Eq.  (4). 


In  Refs.  8  and  9  it  is  shown  that  the  result  for  the  scattering  kernel  as  given  in  Eq.  (4)  is 
adequate  to  describe  existing  experimental  results  provided  appropriate  Ax,  Ay,  and  Az 
where  are  chosen.  The  coordinate  results  from  Eq.  (4)  in  dimensionless  form  are  as  follows: 


P 


(77  ) 


exp 


(6) 


where 


R  = 


X 


S' 

X 


T  =  A  S 

X  XX 


and 


I  (  —  2  R  T  ) 

2T  - - exp  -  (R  +  T  )2 

y  exp  (  -  2  R  T  )  r  ' 

y  y 


[-<VT,)2] 


where 


R  =  yj\2  -  1  S' 

y  v  y 


Ty  =  A  S 
’  y  y 


The  term  Pz  can  also  be  expressed  by  Eq.  (6),  with  the  x’s  replaced  by  z’s. 


2.2  MONTE  CARLO  SOLUTION  OF  P(v',v) 


(7) 


Physically  Eqs.  (6)  and  (7)  give  the  probability  that  a  molecule  with  velocity  components 
vx,  Vy,  and  vz  will  be  reflected  from  a  surface  with  velocity  components  vx,  vy,  and  vz.  Using 
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the  Monte  Carlo  solution  there  are  two  methods  that  can  be  used  to  determine  the  reflected 
velocities  from  Eqs.  (6)  and  (7).  One  of  these  is  the  rejection  technique,  and  the  other  is  the 
direct,  or  inversion,  method.  For  a  brief  description  of  these  methods  see  Ref.  10. 

The  rejection  technique  is  used  frequently  when  involved  mathematical  equations  are 
necessary  for  the  evaluation  of  inverse  probability  functions.  Equation  (6)  and  especially  Eq. 
(7),  which  includes  the  Bessel  function,  IG,  are  time-consuming  to  evaluate.  Since  the  Monte 
Carlo  method  requires  many  samples  to  give  meaningful  results,  the  evaluation  of  any 
involved  numerical  scheme  should  be  avoided.  It  was  therefore  decided  to  use  the  direct 
method  with  a  numerical  approximation  to  the  inverse  scattering  kernel,  as  follows. 

Define  functions  Qx  and  Qy  as 


T 

W  ■  / 

—  oo 

and 

T 

w  -  / 

0 

Since  Px  and  Py  are  normalized,  the  values  of  Qx  and  Qy  vary  from  zero  to  unity.  It  can 
be  shown  that  the  probability  density  of  Qx  and  Qy  is  unity;  in  other  words,  for  a  given  Sx 
and  Sy  any  Qx  or  Qy  between  zero  and  one  is  equally  probable.  This  fact  is  the  basis  of  the 
direct  selection  technique.  By  obtaining  a  random  number  between  zero  and  one  and 
equating  that  to  Q,  one  can  solve  Eq.  (8)  or  (9)  for  T  and  thus  find  the  reflected  velocity 
components  from  the  definitions  of  Tx,  Ty,  Sx,  and  Sy. 

As  stated  above,  a  curve  fit  to  Eqs.  (6)  through  (9)  will  be  used.  Given  Q  and  R,  the 
problem  is  to  find  T  and  therefore  vx,  vy,  and  vz.  It  is  desirable  to  use  a  linear 
approximation,  and  it  was  found  that  by  redefining  variables  a  suitable  form  was  obtained. 
For  example,  R  can  vary  over  all  values,  and  the  curves  are  asymptotic  with  R  as  a 
parameter;  however,  if  one  uses  R2/(0.5  +  R2)  as  a  parameter,  whose  value  is  between  zero 
and  unity,  the  curves  are  more  linear.  The  other  parameter  used,  which  also  has  limited 
variation,  was 

T-R2  /  (0.5  +  R) 

From  the  asymptotic  expansion  of  the  Bessel  function,  IOJ  it  can  be  shown  that  as 
(2  Ry)  *  oo,  Eq.  (7)  reduces  to  Eq.  (6).  This  fact  is  used  to  permit  1)  use  of  the  results  for 


P  r|T 


P  dT 

y  y 


(8) 


(9) 
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Eqs.  (7)  and  (9)  and  2)  inclusion  of  the  results  for  the  x  and  z  components  as  given  by  Eqs. 
(6)  and  (8). 

Based  upon  the  results  of  the  reciprocity  relation  and  the  numerical  curve  fitting 
described  above,  a  computer  subroutine  was  written  to  calculate  the  reflection  from  a 
surface.  This  subroutine,  REFLCT,  is  given  in  Appendix  A. 

2.3  NUMERICAL  RESULTS  FROM  GAS-SURFACE  MODEL 

In  order  to  check  the  reciprocity  gas-surface  model  and  its  implementation  in  REFLCT 
and  also  to  gain  insight  into  the  characteristics  of  reflected  molecular  properties,  a  simple 
parametric  analysis  was  performed.  Assuming  a  molecular  beam  of  argon  with  normal  and 
tangential  velocity  components  of  10,000  ft/sec,  which  is  a  speed  of  14,142  ft/sec  at  an 
incidence  angle  of  45  deg,  the  reflected  properties  were  determined  for  surface  temperatures, 
Ts,  of  300  and  30°K  and  also  for  energy  accommodation  coefficients,  a,  of  0.9  and  0.5. 
Since  multiple  encounters  of  molecules  with  the  surface  are  important  in  the 
accommodation  to  space  chamber  conditions,  each  molecule  was  allowed  to  undergo  10 
collisions  with  the  test  surface.  A  total  of  100,000  molecules  was  followed  for  the  10  surface 
interactions;  thus  1,000,000  gas-surface  encounters  were  calculated  using  the  Monte  Carlo 
technique  for  each  condition. 

In' a  real  enclosure  the  number  of  encounters  a  molecule  will  have  with  surfaces  can  vary 
from  zero  to  an  unlimited  number.  After  a  molecule  enters  the  enclosure  through  an 
entrance  area,  it  is  followed  as  it  makes  encounters  with  various  surfaces  until  it  is  captured 
by  a  surface  (i.e.,  cryopumped)  or  leaves  the  enclosure  through  the  exit  or  entrance.  Also, 
tracking  of  the  molecule  can  be  stopped  after  a  number  of  surface  collisions,  which  can  be 
input  into  the  program.  However,  in  the  present  example  no  surface  capture  takes  place,  and 
no  molecules  can  exit.  The  number  of  surface  encounters  allowed  (10)  is  completely 
arbitrary. 

For  an  isothermal  enclosure,  molecules,  independent  of  their  initial  velocities,  will 
approach,  after  a  sufficient  number  of  surface  collisions,  a  Maxwellian,  fully 
accommodated  velocity  distribution.  Spatially,  a  fully  accommodated  distribution  is 
equivalent  to  a  cosine  distribution. 

The  results  of  this  numerical  study  are  presented  in  Figs.  2  through  7  in  two  different 
ways.  First,  in  Figs.  2  and  3  the  reflected  velocity  distribution  is  shown  for  various  values  of 
surface  collisons,  N,  for  a  =  0.5  and  0.9  and  Ts  =  30  and  300°K.  Notice  that  for  a  =  0.9 
(Fig.  3)  the  velocity  distribution  has  become  a  fully  accommodated  Maxwellian  distribution 
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at  the  surface  temperature  for  10  surface  collisions.  Second,  in  Figs.  4  through  7  the  results 
are  given  as  reflected  spatial  or  angular  distributions.  For  the  case  of  a  -  0.9  and  N  =  10  as 
shown  in  Fig.  7  the  fully  accommodated  cosine  distribution  is  obtained. 


Figure  2.  Reflected  velocity  distribution  for  Ts  =  300°  K  and  a  =  0.5. 


Figure  3.  Effect  of  surface  temperature,  Ts,  on  reflected  velocity 
distribution  for  a  =  0.9. 
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Note:  Radial  coordinate  is  proportional  to  the  number  of 
molecules  reflected  in  a  small  solid  angle. 
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Figure  4.  Reflected  spatial  distribution  for  a  =  0.5  and  Ts  =  300°  K. 
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Note:  Radial  coordinate  is  proportional  to  the 
number  of  molecules  reflected  in  a  small  solid  angle. 
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Figure  4.  Concluded. 


Figure 


AEDC-TR-79-76 


Note:  Radial  coordinate  is  proportional  to 
the  number  of  molecules  reflected  in  a  small 
solid  angle. 
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Figure  7.  Reflected  spatial  distribution  for  a  =  0.9  and  N  =  10. 

One  characteristic  of  the  Monte  Carlo  technique  evident  in  these  calculations  is  the 
random  scatter  of  the  results.  This  is  best  seen  in  Fig.  7,  where  most  of  the  points  are  within 
±5  percent  of  the  mean.  It  can  be  shown  from  statistical  theory  that  the  variation  of  a 
random  process,  such  as  the  Monte  Carlo  method,  about  the  mean  is  a  Gaussian  distribution 
with  a  width  or  mean  error  proportional  to  N_1/2.  Thus,  in  order  to  reduce  the  numerical 
scatter  by  a  factor  of  two,  four  times  as  large  a  sample  would  be  required  (i.e.,  N  =  400,000 
for  the  present  calculation). 

Based  upon  the  results  of  this  limited  analysis  the  following  conclusions  can  be  drawn: 

1.  For  N  =  1  the  reflection  is  close  to  specular,  with  almost  all  reflected  molecules 
having  the  same  tangential  direction  as  the  incident  beam. 

2.  The  reflected  distribution  for  lower  surface  temperature  is  narrower  both  in 
velocity  spread  and  in  spatial  variation  than  that  for  higher  surface  temperature. 
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3.  For  higher  accommodation  coefficients  the  reflection  is  more  diffuse,  and 
equilibrium  conditions  are  approached  with  fewer  surface  collisions  than  for 
lower  values  of  a. 

As  shown  in  Refs.  7,  8,  and  9,  the  reciprocity  gas-surface  scattering  model  is  the  only 
existing  theory  that  agrees  with  available  experimental  data.  In  addition,  the  implementation 
of  the  present  model  in  REFLCT  also  realistically  describes  the  accommodation  of 
molecules  to  the  conditions  of  an  enclosure. 

On  the  basis  of  these  facts,  a  duplicate  of  the  existing  GUERAP  computer  code  was 
modified  to  eliminate  thermal  radiation  peculiar  parts  and  add  the  improved  gas-surface 
model  subroutine  REFLCT.  Changing  the  energy  ray-trace  technique  to  a  molecule  ray  trace 
prompted  changing  the  E  in  GUERAP  to  an  M,  giving  rise  to  the  name  GUMRAP  for  the 
modified  version. 

No  details  on  the  changes  made  to  GUERAP  are  given  here.  However,  a  detailed 
description  of  input  variables  required  to  run  GUMRAP  is  presented. 

3.0  NUMERICAL  INPUTS  TO  GUMRAP 

This  section  contains  a  detailed  definition  of  important  input  variables.  Appendix  B  lists 
all  input  parameters,  a  simplified  description,  and  the  default  values  used  in  the  program  if 
no  overriding  values  are  input.  It  is  suggested  that  one  make  a  sketch  of  the  complete  system 
and  check  all  inputs  before  attempting  to  run  the  program. 

There  are  several  types  of  inputs  to  GUMRAP,  including  the  following: 

1.  Geometric  inputs  of  system  configuration 

2.  Properties  of  surfaces  that  make  up  the  geometric  system 

3.  Properties  of  incident  molecular  beam  in  regard  to  molecular  weight,  velocity, 
direction,  etc. 

4.  Input  data  in  regard  to  random  number  generation 

5.  Options  for  printing  output. 

Most  data  are  input  via  the  namelist  statment,  a  standard  feature  of  FORTRAN. 


17 


A  E  D  C-TR  -79-76 


3.1  GEOMETRIC  INPUTS  (NMGEOM) 

The  first  namelist  is  called  NMGEOM  and  inputs  most  of  the  geometry  of  the  system. 
GUMRAP  can  analyze  systems  of  up  to  30  surfaces  or  boundaries.  Any  system  can  also  be 
divided  into  as  many  as  five  sections.  In  order  to  simplify  the  geometric  input,  seven 
standard-type  surfaces  have  been  defined. 

Each  surface,  identified  by  N,  is  assigned  a  number  and  is  used  to  identify  section 
boundaries,  surface  constraints,  etc.  The  seven  standard  surface  types  will  be  described  first. 
These  seven  surfaces  are: 


Plane 

5. 

Ellipsoid 

Cone 

6. 

Hyperboloid 

Sphere 

7. 

Baffle 

Paraboloid 

The  parameters  for  these  surfaces  as  well  as  optional  standard  constraints  are  input  by 
the  FORTRAN  variable  array  PLTCTR(11,30)  where  up  to  11  values  may  be  input  to 
describe  each  of  the  up  to  30  surfaces.  The  type  of  each  surface  is  identified  by  a  group  of 
arrays  as  follows: 


IPLANE(I)  =  NS  (Surface  number  of  Ith  Plane) 

ICONE(I)  =  NS  (Surface  number  of  Ith  Cone) 

ISHPER(I)  =  NS  (Surface  number  of  Ith  Sphere) 

IPARAB(I)  =  NS  (Surface  number  of  Ith  Paraboloid) 

IELPS(I)  =  NS  (Surface  number  of  Ith  Ellipsoid) 

IHYPB(I)  =  NS  (Surface  number  of  Ith  Hyperboloid) 

IBAFF(I)  =  NS  (Surface  number  of  Ith  Baffle) 

The  inputs  to  PTLCTR  are  detailed  in  Table  I  for  the  seven  standard  surfaces.  Figures  8 
through  16  define  the  various  parameters.  Notice  that  the  baffle  (Figs.  15  and  16)  has  no 
outer  constraint;  special  constraints,  which  will  be  discussed  later,  must  be  used. 
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Table  1.  Definition  of  PTLCTR(J,I\!)  Inputs 


Surface 

Figure 

Number 

PTLCTR(1,N)  = 

Remarks 

Plane  (Disk) 

8 

V  Yc>  Zc>  a,  B,  Y »  R, >  R2 

Unit  Normal  should 
point  into  region. 

Plane  (Rectangle) 

9 

Cone 

10 

Sphere 

11 

Paraboloid 

12 

X1-  V  Z1-  x2-  Y2>  Z2*  V  r2 

Ellipsoid 

13 

Xr  Y1>  Z1-  X2'  Y2>  Z2'  R'  Rc 

Hyperboloid 

14 

Plane  Baffle 

15 

Xc,  Yc>  Zc,  a,  6,  Y,  R,  0. 

Also  Input 

RKNIFE  =  Re 

DBAFF  =  5 

ABAFFD  =  e 

TBAFF  =  t 

Conical  Baffle 

16 

X.,  Yc,  Zc,  a,  6,  Y.  R.  * 

z 


Yc  Zr 


Figure  8.  Plane  (disk)  configuration. 
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Z 

Figure  11.  Sphere  configuration. 


y 


Figure  12.  Paraboloid  configuration. 
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More  general  surfaces  can  be  input  via  the  variable  COEF  where 
COEF  (J,N)  =  Cj  J  =  1,2,  ...  10 

Cj  is  used  in  the  surface  equation,  as  follows: 

CjX2  +  C2Y2  +  C3Z2  +  cxxy  +  c5xz  +  c6yz 

+  c7  X  +  C8Y  +  C9z  =  C10  (10) 

If  needed,  up  to  five  auxiliary  coordinate  systems  can  be  defined  for  referencing  the 
geometric  input  data.  The  particular  coordinate  system  for  surface  N  is  identified  by  the 
variable  ICORD(N)  =  K  where  K  =  0,1... 5  (K  =  0  is  the  system  coordinate  system).  The 
information  relating  each  auxiliary  coordinate  system  to  the  main  system  coordinates  input 
through  the  variable  TRFORM(I,J,K).  For  J  =  1  and  I  =  1,2,3,  TRFORM(I,J,K)  gives  the 
X,Y,  and  Z  direction  cosines  of  the  Kth  system’s  X  axis;  for  J  =  2,  TRFORM(I,J,K)  gives 
the  Y-axis  direction  cosines,  and  for  J  =  3  it  gives  the  Z-axis  cosines.  For  J  =  4, 
TRFORM(I,  J,K)  gives  the  orgin  of  the  Kth  system  as  X,  Y,  and  Z  in  system  coordinates  for  I 
=  1,2,  and  3,  respectively. 

Also  input  through  the  NMGEOM  namelist  is  IDGREE(N).  For  any  surface  input  by 
COEF,  IDGREE(N)  for  that  surface  must  be  specified  as  2,  or  it  will  be  computed  for  the 
basic  surfaces. 


The  last  variable  input  by  means  of  NMGEOM  is  LISTOP(I),  I  =  1,11.  LISTOP  gives 
the  list  of  options,  a  value  of  1,  which  causes  printing,  and  a  value  of  0,  which  suppresses 
that  print  option.  These  options  for  I  =  1,11  are  as  follows: 


LISTOP 


1 .  Summary  of  surface  input  data 

2.  Summary  of  constraint  input  data 

3.  Equations  of  surfaces  and  constraints 

4.  Summary  of  system  configuration 

5.  No  longer  used 

6.  Surface  property  coefficients 

7.  Spatial  maps 
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8.  Radial  map  and  plots 

9.  Spatial/angular  maps 

10.  Table  of  energy  absorbed  by  each  surface 

11.  Combined  cases  computations 


n  =  Unit  Normal  to  Plane  of  Baffle 
(Pointing  Away  from  Cone) 

Figure  15.  Plane  baffle  configuration. 


3.2  CONSTRAINT  INPUTS  (NMCNST) 

Each  surface  input  through  namelist  NMGEOM  must  have  defined  boundaries  or 
constraints.  Several  of  the  standard-type  surfaces  had  standard  constraints  such  as  the  plane 
disk  or  rectangle.  Variables  input  through  namelist  NMCNST  specify  the  type  constraint  or 
define  new  constraints. 
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(Pointing  Away  from  Double  Cone) 

Figure  16.  Conical  baffle  configuration. 

The  variable  NSTRN(N)  is  equated  to  the  number  of  constraints  of  surface  N.  The 
number  of  constraints  per  surface  must  be  four  or  less;  i.e., 

NSTRN(N)  <  4 

Three  type  constraints  can  be  used.  First,  the  standard  as  defined  in  PTLCTR(I,N); 
second,  other  boundary  surfaces  used  as  defined  in  namelist  NMGEOM;  and  third,  a  new 
constraint  surface  defined  by  the  parameter  PTLCSP(I,M).  A  plane  or  conical  special 
constraint  surface  can  be  input  by  means  of  PTLCSP(I,M)  by  the  same  format  as  used  in 
PTLCTR(I,N)  and  defined  in  Table  1.  However,  no  constraints  are  required  for 
PTLCSP(I,M)  since  it  itself  is  a  constraint.  The  first  special  constraint  is  identified  by  M  = 
1,  the  second  by  M  —  2,  etc.  There  is  no  necessary  relation  between  the  surface  number  and 
constraint  number. 
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The  special  constraint  surfaces  can  be  input  using  the  auxiliary  coordinate  systems  that 
were  discussed  and  could  be  introduced  in  NMGEOM.  If  auxiliary  systems  are  used,  they 
must  be  identified  by  the  variable  ICORSP(M).  This  variable  must  be  equated  to  the  number 
of  the  particular  auxiliary  system  used  to  input  the  special  constraint. 

The  variable  KSIDSP(I,N)  defines  the  l'1'  constraint  of  surface  N  as  to  type  as  follows: 

1.  For  an  existing  surface  defined  in  NMGEOM,  set  KSIDSP  equal  to  the 
constraint  surface  number. 

2.  For  the  standard  constraint  as  input  by  PTLCTR,  set  KSIDSP  equal  to  zero. 

3.  For  a  special  constraint  defined  by  PTFCSP(I,M),  set  KSIDSP  equal  -M. 

One  must  also  specify  on  which  side  of  a  constraint  the  surface  of  interest  lies.  This 
information  is  input  via  the  variable  KSTRT(J,N)  for  the  Jth  constraint  on  surface  N,  as 
defined  in  Table  2. 


Table  2.  Definition  of  KSTRT(J,!\!) 


Type  Constraint 

KSTRT(J.N)  = 

Plane 

1 ,  when  surface  normal  of  constraint  points 
away  from  surface 

2,  when  surface  normal  points  toward  surface 

Cone ,  Sphere 

Paraboloid, 

Ellipsoid 

3,  when  surface  is  on  concave  side  of  constraint 

4,  when  surface  is  on  convex  side  of  constraint 

Hyperboloid 

3,  when  surface  is  on  convex  side  of  constraint 

4,  when  surface  is  on  concave  side  of  constraint 

General  Surface 

F (X, Y, Z)  =  0 
where  „  2  2 

F  =  +  C3Z 

+  C4XY  +  C5XZ  +  CfeYZ 

+  c7x  +  c8Y  +  c9z  -  c10 

3,  when  F  <  0  on  surface  of  interest 

4,  when  F  >  0  on  surface  of  interest 

Each  special  or  standard  constraint  requires  computer  memory  in  order  to  input  the 
constraint  surface.  Each  constraint  defined  in  this  way  lowers  by  one  the  number  of  surfaces 
which  can  be  used  to  define  the  system  boundary.  For  example,  a  simple  box  enclosure, 
consisting  of  six  plane  surfaces,  if  input  via  the  standard  rectangle  input  with  standard 
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constraints,  would  require  four  constraints  on  each  face  as  well  as  the  six  faces.  This  total 
would  require  30  surfaces,  which  is  the  maximum  the  present  version  of  GUMRAP  will 
handle.  This  means  that  each  surface  is  used  five  times,  once  as  a  surface  and  four  times  as 
constraints;  this  is  unnecessarily  redundant.  If  each  surface  were  defined  and  then  used  as  a 
constraint,  only  six  out  of  thirty  memory  locations  would  be  required.  Thus,  the  way  in 
which  the  constraints  are  input  becomes  very  important  in  utilizing  the  full  capability  of  the 
program. 

3.3  SURFACE  PROPERTY  INPUTS  (NMCOAT) 


Up  to  25  different  type  surfaces  can  be  defined.  The  definition  of  each  type  surface  and 
the  identification  of  each  surface  type  is  input  through  namelist  NMCOAT. 

The  surface  properties  are  input  by  use  of  the  variable  COAT(I,J),  which  is  the  Ith 
property  of  the  Jth-type  surface.  The  elements  of  COAT  are  defined  as  follows: 


COAT(l,J) 

COAT(2,J) 

COAT(3,J) 

COAT(4,J) 

COAT(5,J) 

COAT(6,J) 


the  probability  that  a  molecule  will  be  reflected  from  surface  type  J 
according  to  the  reciprocity  model; 

the  probability  that  a  molecule  will  be  reflected  from  surface  type  J 
specularly; 

the  probability  that  a  molecule  will  be  reflected  diffusely  with  no 
change  in  speed; 

surface  temperature  (°K); 

the  tangential  energy  accommodation,  ax  and  ay; 
the  normal  energy  accommodation,  ay. 


The  last  three  variables  can  be  omitted  if  COAT(l,J)  =  0. 


The  variable  ICOAT(N)  identifies  each  surface  N  with  a  particular  type  surface  J;  i.e. , 

ICOAT(N)  =  J 


3.4  HEADING  AND  ISEE1) 

Any  number  of  alphanumeric  cards  may  be  input  after  the  end  of  the  NMCOAT 
namelist.  These  cards  will  be  printed  out  twice  as  headings.  The  heading  cards  must  be 


27 


AEDC-TR-79-76 


terminated  by  the  word  EXTERNAL  (starting  in  column  1).  Next  an  odd  integer  called 
ISEED  is  input  in  columns  1 1  through  20  subject  to 

ISEED  <  231 

ISEED  is  used  to  start  the  pseudorandom  number  sequence  generation  subroutine  used 
by  the  Monte  Carlo  method.  For  the  same  value  of  ISEED  the  sequence  of  random  numbers 
produced  will  be  repeated. 

3.5  INPUT  DATA  DEFINING  SYSTEM  SECTIONS  (NMSYST) 

As  mentioned  earlier,  a  system  can  be  divided  into  as  many  as  five  sections.  The 
definition  of  these  sections  is  input  through  the  namelist  NMSYST.  As  stated  previously, 
system  boundary  surfaces  are  input  through  namelist  NMGEOM,  and  the  constraints  are 
input  through  namelist  NMCNST;  in  addition,  it  is  necessary  to  input  data  to  describe  on 
which  side  of  a  boundary  surface  a  region  or  section  lies.  This  is  accomplished  in  this  group 
of  inputs. 

The  surfaces  that  compose  a  given  section  are  identified  by  the  variable 

ISFACE(J,I)  =  N 

where  ISFACE(J,I)  describes  which  system  surface,  N,  defines  the  Jth  surface  of  section  I. 

The  variable  KCOEF(J,I)  indicates  which  side  of  the  Jth  boundary  surface  of  section  I 
the  section  is  on.  KCOEF  can  have  values  of  0  or  1,  as  shown  in  Table  3. 

When  a  system  is  partitioned  into  smaller  regions  or  sections,  the  common  surface 
between  adjacent  sections  will  have  been  input  as  a  boundary  surface.  The  variable 
ISFTP(I,  J)  indicates  how  the  various  sections  are  interconnected.  This  variable  is  defined  as 

ISFTP(I,J)  =  AJ 

where  AJ  is  the  change  in  section  number  in  going  through  the  Ith  surface  of  section  J  into  an 
adjacent  section,  J  +  AJ.  Notice  that  AJ  may  be  positive  or  negative.  The  default  value  is  0, 
which  indicates  a  reflection  back  into  the  original  section  (surface  1  is  not  an  exit). 

3.6  INPUT  DATA  DEFINING  MOLECULAR  SOURCE  (NMEXTL) 

The  last  namelist  (NMEXTL)  is  used  to  input  data  concerning  the  molecular  source, 
entrance  and  exit  apertures,  and  miscellaneous  control  parameters. 
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Table  3.  Definition  of  KCOEF(J,l) 


Type  Surface 

KCOEF ( J , I) 

Plane 

0,  when  normal  of  section  surface  J  points 
away  from  section  I 

1,  when  normal  points  toward  section. 

Cone,  Sphere, 
Paraboloid, 

Ellipsoid 

0,  when  section  I  is  on  concave  side  of  section 
surface  J 

1,  when  section  is  on  convex  side  of  surface 

Hyperboloid 

0,  when  section  I  is  on  convex  side  of  section 
surface  J 

1 ,  when  section  is  on  concave  side  of  surface 

General  Surface 
F(X,Y,Z)  =  0 

0,  when  F  <  0  in  section  I 

1 ,  when  F  >  0  in  section  I 

In  order  for  the  Monte  Carlo  ray  trace  technique  to  be  valid,  molecules  must  be 
introduced,  one  at  a  time,  in  a  realistic  fashion  into  the  system.  One  of  the  surfaces 
composing  the  geometric  systems  can  be  designated  the  entrance  by  the  variable  ISENT 
equated  to  the  surface  number  of  the  entrance  aperture.  The  test  molecules  enter  uniformly 
over  the  entrance  aperture  with  a  direction  and  divergence  given  by  parameters  described 
below  and  shown  in  Fig.  17. 

Figure  17  depicts  a  typical  divergent  molecular  beam  passing  through  a  circular  entrance 
aperture.  The  aperture  shown  is  placed  on  the  X-Y  plane  of  the  system  coordinates, 
although  this  is  not  necessary.  The  molecule  enters  the  aperture  by  the  randomized  Monte 
Carlo  technique.  The  entrance  aperture  can  be  either  circular  or  rectangular  . 

The  direction  of  the  molecular  beam  as  a  whole  is  defined  by  the  elevation  angle,  0, 
which  is  input  via  the  variable  SUNAGD  and  the  azimuthal  angle,  9,  which  is  input  via  the 
variable  SUNAZD.  The  divergence  of  the  beam  (half  angle  6)  is  input  via  the  variable 
SUNDVD.  All  angles  are  in  degrees. 

The  elevation  angle,  <t>,  is  defined  as  the  angle  between  the  beam  axis  and  the  Z  axis. 
Notice  from  Fig.  17  that  if  vz  <  0  then  (/>  <  90  deg.  The  aximuthal  angle  9  is  the  angle 
between  the  projection  of  the  beam  axis  on  the  X-Y  plane  and  the  X  axis.  Other  parameters 
that  are  input  through  the  namelist  NMEXTL  are  given  in  Appendix  B. 
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Figure  17.  Entrance  aperture  and  molecular  beam  geometry. 

3.7  SAMPLE  INPUT  DECK 

A  listing  of  the  card  deck  necessary  to  run  GUMRAP  on  an  IBM  370/165  computer  is 
shown  in  Appendix  C.  Notice  that  for  the  namelist  input  format  only  the  first  element  in  an 
array  name  need  be  given;  for  example: 

PTLCTR(1,8)  =  2*. 5,  2.222,  0.,etc. 


is  equivalent  to 


PTLCTR(1,8)  =  .5 
PTLCTR(2,8)  =  .5 
PTLCTR(3,8)  =  2.222 
PTLCTR(4,8)  =  0.0 
etc. 
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Each  namelist  starts  with  an  ampersand  (&)  in  colunm  2  followed  immediately  by  the 
namelist  name  and  then  a  blank;  after  this,  data  items  follow  on  the  same  or  additional 
cards,  separated  by  commas.  Subsequent  data  cards  must  start  with  a  blank  in  column  1, 
with  data  starting  in  column  2.  A  namelist  is  terminated  by  an  “ampersand  end”  (&END), 
also  starting  in  column  2.  The  order  of  data  in  a  given  namelist  is  immaterial.  However,  the 
order  of  the  namelists  must  be  followed  as  given. 

4.0  SUMMARY 

The  General  Unwanted  Energy  Rejection  Analysis  Program  (GUERAP)  was  used  as  a 
basis  for  a  molecular  flow  analysis  program  (GUMRAP).  The  system  geometry  and  Monte 
Carlo  ray  trace  portion  of  GUERAP,  which  was  the  greater  part,  was  combined  with  a  new 
Monte  Carlo  gas-surface  interaction  program,  REFLCT. 

It  was  shown  that  REFLCT  can  be  used  to  calculate  the  resulting  conditions  for  a  highly 
nonequilibrium  molecular  beam  entering  an  enclosure,  and  that  the  Maxwellian  velocity 
distribution  as  well  as  the  fully  accommodated  cosine  distribution  is  approached  after  10 
collisions  with  the  chamber  for  an  accommodation  coefficient  of  0.9. 

All  input  conditions  were  defined  so  that  GUMRAP  can  be  run  on  an  IBM  370/165 
Computer  for  complex  internal  flow  geometries. 
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APPENDIX  A 

LISTING  OF  SUBROUTINE  REFLCT 


SUBROUTINE  REFLCT ( V»N) 

GIVEN  A  MOLECULE  WITH  VELOCITY  VECTOR  (V)  INCIDENT  UPON  A  SURFACE,  N» 
THIS  SUBROUTINE  RETURNS  IN  V  THE  REFLECTED  VELOCITY  VECTOR  BASED 
UPON  THE  PECIPROCIRY  RELATION,  VELOCITY  IN  FT/SEC. 

I  SEED  in  COMMON/SETl/  is  USED  IN  FUNCTION  RN2  TO  GENERATE  A  RANDOM  Q, 
WMOLT  IN  COMMON/SE T2/  IS  THE  GAS  MOLECULAR  WEIGHT. 

COAT (7, N)  AND  COAT <8, N)  IN  COMMON/SE T3/  ARE  AX  AND  AY  RESPECTIVELY. 

COAT  (4, V)  IN  COMMON/SE T 0/  IS  THE  SURFACE  TEMPERATURE  IN  DEGREE  K  ,TS, 
OF  SURFACE  N. 

ugasc  is  the  universal  gas  constant  in  (ft#*2)/(sec**2> ioeg  k> 

DIMENSION  F  T  ( 2  7 , 1 1 ) , FIT (27,5) » V  ( 3 ) » A  ( 3 ) 

COMMON  /SET  I /  I  SEED, KEMI T  *KSID,KSTYPE ,NSEC , KSEC , IERR , NDI aG» I POLE « 

ikssec,npole»ksecd»irun 

COMMON  /SET2/  KS ID a, COEF (10,30) ,KSTRT (4 , 30 ) , CSTRT ( 1 0 ,4 ) ,NQUAR, 

1 PTLCTR (11,30) , NSTRT (30) ,1DGREE(30> ,KSIDSP (4, 30) ,ISTRT(4,30) , 

2LSURF  C,ICORD(30)  , TRFORM ( 3,4,5) , LSEC , LSSEC , TCONO , VEL 1NT , WMOLT 
C0MM0N/SET3/  RAY (b) ,SURN0R(3) , COAT (4 1,25) , E , PH  I , THE , RA YR (6) , 
1ESAVE,ETHERM(30) , I EMOOL  <  25 ) ,LlSTOP(il) ,ICOAT(30) 

DATA  IjGASC/89523  « / 

EQUIVALENCE  < F I T ( 1 , 1 ) , F T ( i , 7 ) ) 

DATA  FT/ 

2  ,0260,  .1427,  .2023,  ,2489,  .2888,  .3246,  .4032,  .4724,  .5364, 

3  .5972,  .6563,  .7147,  .7132,  ,8326,  ,8936,  .9572,1.0246,1.0973, 

4 1. 1 7 74, 1.2686, 1. 3 ?74, 1.51 74, 1,5892, I, 6773, 1,794 1, 1,9779, 2. 4930, 

5-. 0400,  ,0713,  .1325,  .1804,  .2214,  .2582,  .3390,  .4101,  .4759, 

6  .5384,  .5992,  .6692,  .7193,  .7803,  .8430,  .9083,  .9775,1.0521, 

71.1344,1.2280,1.3396,1.4833,1,5570,1.6473,1.7669,1.9552,2.4823, 

8". 1030,  ,0056,  ,0b88,  .1184,  .1609,  .1990,  .2825,  .3562,  ,4242, 

9  .4888,  .5516,  .6136,  .b756,  .7386,  .8032,  .8706,  .9420,1.0188, 

A  1.1 035, 1.  1 998, 1.3 145, 1,46 19, 1,53 75, 1.6300, 1.7525, 1.9449, 2, 4821, 

H-. 1840, -.0637,  .0023,  ,0542,  .0987,  .1385,  .2257,  .3025,  .3734, 

C  .4407,  .5060,  .5705,  ,6349,  .7002,  .7673,  .8371,  .9109,  .9903, 

01.0777, 1 , 1769,1,2949, 1 .4462, 1 .5236, 1 .6183,1 . 7435 , 1 . 9397 ,2.4850 , 

E-. 2670, -.1412, “.0/08, -.0158,  .0313,  .0734,  .1655,  .2466,  .3213, 

F  ,3921,  ,4607,  .5282,  ,5956,  ,6638,  ,7337,  .8064,  .8830,  .9652, 

Gl. 0556, 1.1579, 1.2 792, 1.4344, 1.51 36, 1.61 03, 1.7379, 1.9374, 2. 4895, 

H-. 3690, -.2317, -.1551 ,-.0956, -.0447,  ,0007,  .1000,  .1869,  .2667, 

I  .3421,  .4149,  .4863,  .5574,  ,6291,  .7023,  .7782?  .8579,  .9432, 

J 1.0367, 1 . 1422,1,2660,1.4256,1.5064,1.6050,1.7347,1.9371,2.4948/ 

DATA  FIT/ 

K-. 4890, -.3424, -.2564, -,189b, -.1327, -.0821,  .0276,  .1229,  ,2097, 

L  .2912,  .3694,  .4456,  .5211,  .5967,  .6737,  .7530,  .8360,  .9244, 

Ml .0209, l .1294,1.2671, 1.4193,1  .5015,1.6016,1.7333,1.9381,2.5006, 
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N-.  6740  4854, -.  38 1 5, -.  3020, -.  23'S  1  1762,-.  0S05,  .0567,  .1528? 

0  ,2419,  .3263,  .4079,  .4881,  .5679,  .6485,  .7312,  .8173,  .9087? 
PI. 0080, 1.1 192, 1.2496, 1,4 147, 1.4982, 1.5997, 1.7330 *1.9400, 2. 5069, 
Q-.9430,-. 6725, -.53 12, -.4284, -.3449, -.2734,-. 12S9.-.0047,  .1015, 

R  .1982,  .2888,  .3/54,  .4597,  ,5432,  ,6271,  .7127,  .8015,  .8954, 

S  .9971,1.1107,1.2435,1.4112,1.4960,1.5988,1.7336,1.9427,2.5140, 
r-1 .3420, -.8577, -.06 19, -.5330, -.4335,-, 3510,-. 1861 ,-,0542,. 0895, 

U  .1618,  .2570,  .3*474,  .4351,  .5215,  .6080,  .6961,  ,7 8/2,  .8833, 

V  .9872,1.1031,1,2382,1.4085,1.4944,1.5986,1.7350,1.9463,2.5225, 
W- 1.53 70, -.9522, -.^379, -.5994, -.4935, -.4062, -.2329, -.0951,  ,0231, 
X  .1292,  .2275.  .3209,  .4ili,  .5000,  .5889,  ,6791,  .7725,  .8708, 

Y  .9769,1.0951,1.2329.1.4062,1.4935,1.5994,1.7379,1.9522,2.5355/ 

=. WMOLT 
FS  =  COAT (4,N) 

\/U=SQRT (2.*UGASC*rs/WM) 

A(l)  =  COAT (  7  ,  N ) 

A  ( 2 )  =  COAT ( 8, N ) 

A ( 3 )  =  COAT (7,N) 

00  4  K  — 1,3 
SSV (K) /VW 

R  =  SORT ( A (K ) **2  -  1.0)#S 
CALL  RN2  ( I  SEEl) , Q ) 

FR  =  1 . 

IF  (K.EQ.2)  FR=R*R/(.5  +  R.*R) 

IF  (CJ.LE.0.1)  GO  TO  1 

IF  (Q.GT.0.9)  GO  TO  2 

IX=20*Q 

I=IX  +4 

X320 . *0- 1 X 

GO  TO  3 

1  1X=50#Q 
1=1X41 
X=50.*0-IX 
GO  TO  3 

2  IX=50*0 
I=IX-23 
X=50.*Q-IX 

3  IF (FR.EQ.1.0)  GO  TO  5 
I  Y=1 0»FR 

J=IY  +1 
Y=10.#FR-I Y 

FTO=  ( 1  .  -Y )  *  (  U.-X)*FTd,J)4X*FT<l4l,J)  ) 
t* Y*( <1.-X)*FT (I *  J  +  1 ) >X*FT  <1*1, J+l)  ) 

T  =  FTO  ♦  <R*R/<0.5  ♦  R )  ) 

GO  TO  6 

5  FT0=  (l.-X)*FTd,ll)4X*FT(l4l,H) 

T=FTO-.S4« 

6  V (K)  =  T#VW/A ( K ) 

4  CONTINUE 
RETURN 
ENO 
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APPENDIX  B 

LISTING,  DEFINITION,  AND  DEFAULT  VALUES  OF  INPUT  VARIABLES 
NAMELIST  —  NMEGEOM 


Variable 

Definition 

Default  Value 

IPLANE(I) 

Surface  number  of  the  Ith  plane. 

0 

ICONE(I) 

Surface  number  of  the  Fh  cone. 

0 

ISPHER(I) 

Surface  number  of  the  Fh  sphere. 

0 

IPARAB(I) 

Surface  number  of  the  Fh  paraboloid. 

0 

IELPS(I) 

Surface  number  of  the  Ilh  ellipsoid. 

0 

IHYPB(I) 

Surface  number  of  the  Ith  hyperboloid. 

0 

IBAFF(I) 

Surface  number  of  the  Fh  baffle. 

0 

PTLCTR(J,N) 

Jth  parameter  describing  standard  surface  N  (see 
Table  1). 

0 

COEF(J,N) 

Coefficient,  Cj,  for  special  surface  N  defined  by 
Eq.  (32). 

0 

RKNIFE 

Re  (see  Figs.  15  and  16). 

0 

TBAFF 

t,  thickness  of  baffle  (see  Figs.  15  and  16). 

0 

DBAFF 

<5  (see  Fig.  15). 

0 

ABAFFD 

e,  baffle  blade  angle  (see  Figs.  15  and  16). 

0 

IDGREE(N) 

Degree  of  equation  for  surface  N  =  2  for  special 
surface. 

Computed  for 
standard  surfaces 

TRFORM(I,J,K) 

Rotation  matrix  and  origin  of  auxiliary 
coordinate  system  K  (see  Section  3.1). 

ICORD(N) 

=  K,  coordinate  system  used  for  input  data  of 
surface  N. 

0  (System 
coordinates) 

LISTOP(I) 

=  1  for  any  option  wanted  (see  Section  3.1). 

0  (No  output) 
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NAMELIST  —  NMCNST 


Variable 

Definition 

Default  Value 

NSTRT(N) 

Number  of  constraints  on  surface  (N),  <  4, 
must  be  entered  for  all  surfaces. 

0 

KSIDSP(I,N) 

Definition  of  the  I,h  constraint  type  for  surface 

N  (see  Section  3.2). 

0 (Standard 
defined  by 
PTLCTR) 

KSTRT(I.N) 

Type  of  the  Ith  constraint  of  surface  N  (see 
Table  2). 

0 

PTLCSP(I.N) 

Special  Fh  constraint  for  surface  N  (plane  or 
cone)  by  same  format  as  PTLCTR  (Table  1). 

0 

ICORSP(M) 

=  K,  coordinate  system  used  for  input  data  of 
special  constraint,  M. 

0  (System 
coordinates) 

NAMELIST  - 

-  NMCOAT 

Variable 

Definition 

Default  Value 

COAT(U) 

Fh  properties  of  surface  type  J  (see  Section  3.3). 

0 

ICOAT(N) 

=  J,  surface  type  of  surface  N. 

0 

After  NMCOAT,  cards  containing  a  title  to  be  printed  on  the  output  can  be 
included.  Next  the  word  EXTERNAL  starting  in  column  1  must  be  included, 
followed  by  the  value  of  ISEED  (See  Section  3.4). 

NAMELIST  - 

-  NMSYST 

Variable 

Definition 

Default  Value 

ISFACE(J,I) 

=  N,  the  surface  number  of  the  Jth  surface  of 
Section  I. 

0 

KCOEF(J,I) 

coefficient  indicating  which  side  of  the  Jth 
surface  section  I  is  on  (see  Table  3). 

0 

ISFTP(J,I) 

=  AJ,  change  in  section  number  for  molecules 
exiting  through  surface  I  of  Section  J;  may  be 
+  or  -. 

0 
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NAMELIST  - 
Variable 

NSEC 

LEMIT 

LRFN 

ISEEDP 

NPRINT 

RMAP 

KPRINT 

ISEXT 

JPRINT 

SUNAGD 

SUNAZD 

SUNDVD 

I  SENT 


NMESTL 


Definition 


Default  Value 


Number  of  sections  in  the  system  <  5.  1 

Number  of  molecules  that  enter  the  system.  0 

Limit  on  number  of  molecular  surface  20 

encounters. 

After  every  ISEEDP  molecules  enter  the  system,  100 


the  value  of  ISEED  is  printed  out;  these  ISEED 
values  can  be  used  to  continue  or  repeat  the 
computation. 

Limit  on  number  of  molecules  for  which  the  ray  200 

tracing  will  be  printed  in  arriving  at  the  exit 

aperture. 

Radius  of  the  molecular  distribution  map  at  the  0 

exit  aperture. 

Input  for  printing  ray  trace  history  of  rays  that  0 

did  not  reach  the  exit  aperture.  Set  any  positive 
integer  and  print  will  occur  after  first  KPRINT 
molecules.  If  =0  no  such  printout  will  occur. 

Exit  aperture  surface  number. 


Number  of  molecules  with  one  line  summary  1 

printout. 

Elevation  angle  of  incident  beam,  <j>  (degrees);  0 

see  Fig.  17. 

Azimuth  angle  of  incident  beam,  d  (degrees);  see  0 

Fig.  17. 

Divergence  angle  of  incident  beam,  <5  (degrees);  0 

see  Fig.  17. 

Entrance  aperture  surface  number. 
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NAMELIST  —  NMESTL  (Concluded) 


Variable 

Definition 

Default  Value 

IENTTP 

=  0  for  circular  entrance, 

=  1  for  rectangular  entrance. 

0 

IESTTP 

=  0  for  circular  exit, 

=  1  for  rectangular  exit. 

0 

KSECST 

Section  number  in  which  molecules  originate. 

1 

ICSENT 

Coordinate  system  for  the  entrance  section. 

0 

VELINT 

Initial  velocity  of  molecules  entering  system 
(ft/sec). 

1000 

TCOND 

Condensation  temperature  of  molecules  (°K). 

77. 

WMOLT 

Molecular  weight  of  molecules. 

28. 
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APPENDIX  C 

SAMPLE  GUMRAP  INPUT  DECK  FOR  IBM  370/165 

/^PRIORITY  5 

//ASH05919  J08  ( APO  * 

//  VPV0020J*()l  ,VJ2K- 15A)  ,02489KINSLOW,rtSGLEVEL=l  , 

//  CLASS5  *  *  T1ME=(  *  1  b  )  *PEGl0U=190K 

//  EXEC  F0RTEPDS,PGON0  =  \/9V00203 

//GO  »  F  T  OSH  001  DO  UNI  r  =  .-JOPK,  5PACt=(CYLf  (9,1)  ,  PLSE  »  ,  POUND )  , 

//  DCB= ( P£CFM5VGS  *LPFCL  =  84  »8LKS I ZE=844 ) 

//GO • F 1 04F  001  00  * 

MMGEOM 

IPLANE<1)=1,3*11, 

PTLCTP (1 , 1 ) =.5, .5,0. ,0. ,0.  ,  1 . , 1  .  ,  l  .  ,  1  .  ,0.  ,0.  , 

COEF ( 1 ,2) =7*0. , 1 .  ,-l .  ,  i .  , 

COEF<l»3)=7*0.,l.,-l.,u., 

COEF (1»4)=6#0.,1.,0«*0.,0.» 

COEF (l»b)=6*0»,l«,0*»0.,l#» 

COEF (1 ,6) =7*0. , 1. , ] . ,3.222, 

COEF ( 1 ,7>=7#0. , 1 . , 1 . ,2.222 » 

PTLCTP ( 1 »  8)  =  .5, .5,?. 222,0. , 0 . , - l . , 1 . , 1 . ,  1  . , 0 .  , 0 .  , 

COEF  (l,9)=61*0e,l<>,0»,0«,0»5 
COEF ( 1 , 10) =6*0 . , 1 .  .0. »u.  ,  1  .  , 

PTLCTP  (1 , 11 ) =.5,1 ,61U ,1 . 11 l ,0. .0. , l .  ,1  .  ,1.  ,  1  .  ,0.  ,0.  , 

I0GREE ( 1 )=1 ,6*2,1 ,2,2,  i  , 

LIST OP(I)=l, 1,1, 1,1, 1,1, 1,1, 1,0, 

S.ENO 

&NMCNST 

NSTWT ( 1 ) =1 1 *4, 

KSIDSP ( 1 , 1 ) =2, 3,4,5, 

KSIOSPl 1 ,2) =1 ,4,5,6, 

KSIOSPU  ,31  =  1,4,5,  7* 

KSIOSP (1,4)31,2,3,11, 

KSI0SP(1,5)=1,2,3,11, 

KSIOSPl  1 ,6) =8,9, 10,2, 

KS I DSP (l,7)=d,9,10,3, 

KSIOSP ( 1 ,8) =6, 7,9, 10, 

KSIOSP (1,9) =6,7,8, 1  1, 

KS I OSP ( 1 , 10) =6, 7,8, 1 1  , 

KS I  DSP (1,11)=2*3?4,5, 

KSTPT ( 1 , 1 )  =  1 ,2,2, 1  , 

KSTPT ( 1 ,2) =2,2, 1 , 1 , 

KS  TPT ( 1 ,3) =2,2, 1  ,  1  , 

KSTPT ( 1 ,4) =2, 1 ,2, 1  , 

KSTPT ( 1 ,5) =2, 1 ,2, 1  , 

KSTPT ( 1 ,6) =2,2, 1 , 1 , 

KSTPT (1,7) =2, 2, 1,1, 

KSTPT  ( 1  ,.8)=1 ,2,2,  l , 

KSTPT (1,91=1,2,2,2, 

KSTKT<1,10)=1,2,2*2, 

KST»T (1,111=1,2,2,1, 

S.ENO 
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&NMC04T 

COA  r«l*l)=1.0,0.»0.»j00.,0.5*0.b» 

CO  A  T(l*2)=0.»l«0»0.0» 

ICOAT  U) =2,1 , l»2»2, 1  *  1 *2*2*2»2» 

&ENO 

TEST  CASE  FOR  1X1X1  CHEVRON  BAFFLE  DESIGN 
PROJECT  MO.  V32K-1SA»GUMRAP  MOLECULAR  ANALYSIS 
KINSLOW,  AUGUST  197R. 

EXTERNAL 

13131313 

&NMSYST 

ISEACEU  .1)=1*2»3»4*S»11» 

I SF  ACE  ( 1  »,2 )  =  1 1  *  6  *  7 » B  *  R 1 1 1.)  t 
KCOEF (l.l)=lt0»l»l»0*J* 

KCOEF (l»?)=l»0.1.1.i»Ut 
ISFTP(1. 11=5*0,1* 

ISFTP ( 1 ,2) =-l » 

S.END 

5. NMEXTL 
NSEC=2» 

LEM  I  T  =  5(J  * 

RMAP=.b, 

LRFN=aO, 

KPRINT  =  1 , 

ISEAT=B, 

JPRINT= 1 U  * 

NPR IN  T =50  » 

SUNAGO=iBO. * 

SUNDVU=90. , 

VELINT=20000.  ♦ 

WMOLT=40.  * 

ISENT=1 » 

IENTTP=1* 

IEXTTP=1 , 

6. ENO 

//GO.FT04F 002  DO  * 

/■» 

/# 
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Ax,  Ay,  etc. 

N 

-> 

n 

P(v',v) 

Px(vx>vx)> 

P y( Vy > Vy) ,  etC. 

R 

Rj,R2 

Re 


Rx 

Ry 

S,SX,  etc. 
Ts 
Tx 
Ty 

t 

* 

V 


V  x  ,  V  y  5  etc. 

X,  Y,  Z 
a 


ax,  ay,  etc. 


NOMENCLATURE 

Constant  >  related  to  ax,  ay  by  Eq.  (5) 

Number  of  surface  collisions  in  Section  2.3,  surface  number  in  Section  3 
Unit  normal  vector  with  direction  cosines  (a,(S,y) 

Scattering  kernel 

Scattering  kernel  in  corresponding  coordinate  direction 

Gas  constant  per  unit  mass,  radius  in  Figs.  13,  15,  16 
R.adii 

See  Fig.  16 

VaFTs' 

Va^-  i  s; 

Dimensionless  velocity,  vA/2RTs,  vx/V2RTs,  etc. 

Wall  or  surface  temperature 
AXSX 

AySy 

See  Fig.  15 

Reflected  molecular  velocity 
Incident  molecular  velocity 
Velocity  component  in  corresponding  coordinate 
Coordinate  axis  (see  Fig.  1) 

Thermal  accommodation  coefficient 

Accomodation  coefficient  in  corresponding  coordinate  direction 
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See  Fig.  15  or  Fig.  17 
See  Fig.  15  or  Fig.  16 
See  Fig.  1  or  Fig.  17 
See  Fig.  16 
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